arXiv: 1504.05162v2 [hep-ex] 26 Oct 2015 


EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 



A 

T 

L 

A 

S 



CERN-PH-EP-2015-065 

Submitted to: Physical Review D 


Search for massive, long-lived particles using multitrack 
displaced vertices or displaced lepton pairs in pp collisions at 
y/s = 8 TeV with the ATLAS detector 


The ATLAS Collaboration 


Abstract 


Many extensions of the Standard Model posit the existence of heavy particles with long lifetimes. 
This article presents the results of a search for events containing at least one long-lived particle that 
decays at a significant distance from its production point into two leptons or into five or more charged 
particles. This analysis uses a data sample of proton-proton collisions at y/s = 8 TeV corresponding 
to an integrated luminosity of 20.3 fb' 1 collected in 2012 by the ATLAS detector operating at the Large 
Hadron Collider. No events are observed in any of the signal regions, and limits are set on model 
parameters within supersymmetric scenarios involving R -parity violation, split supersymmetry, and 
gauge mediation. In some of the search channels, the trigger and search strategy are based only 
on the decay products of individual long-lived particles, irrespective of the rest of the event. In these 
cases, the provided limits can easily be reinterpreted in different scenarios. 


© 2015 CERN for the benefit of the ATLAS Collaboration. 

Reproduction of this article or parts of it is allowed as specified in the CC-BY-3.0 license. 




Search for massive, long-lived particles using multitrack displaced vertices or 
displaced lepton pairs in pp collisions at yjs = 8 TeV with the ATLAS detector 

The ATLAS Collaboration 

Many extensions of the Standard Model posit the existence of heavy particles with long lifetimes. 

This article presents the results of a search for events containing at least one long-lived particle 
that decays at a significant distance from its production point into two leptons or into five or more 
charged particles. This analysis uses a data sample of proton-proton collisions at y/s = 8 TeV 
corresponding to an integrated luminosity of 20.3 fb _1 collected in 2012 by the ATLAS detector 
operating at the Large Hadron Collider. No events are observed in any of the signal regions, and 
limits are set on model parameters within supersymmetric scenarios involving A-parity violation, 
split supersymmetry, and gauge mediation. In some of the search channels, the trigger and search 
strategy are based only on the decay products of individual long-lived particles, irrespective of the 
rest of the event. In these cases, the provided limits can easily be reinterpreted in different scenarios. 

PACS numbers: 12.60.Jv, 13.85.Rm, 14.80.Ly, 14.80.Pq 


I. INTRODUCTION 

Several extensions to the Standard Model (SM) pre¬ 
dict the production at the Large Hadron Collider (LHC) 
of heavy particles with lifetimes of order picoseconds to 
nanoseconds (e.g., see Ref. [1] and references therein). 
At the LHC experiments, the decay of a long-lived par¬ 
ticle (LLP) with lifetime in this range could be ob¬ 
served as a displaced vertex (DV), with daughter par¬ 
ticles produced at a significant distance from the inter¬ 
action point (IP) of the incoming proton beams. Parti¬ 
cle decays may be suppressed by small couplings or high 
mass scales, thus resulting in long lifetimes. An exam¬ 
ple of a small-coupling scenario is supersymmetry with 
A-parity violation (RPV) [2, 3]. The present (largely 
indirect) constraints on RPV couplings allow the decay 
of the lightest supersymmetric particle (LSP) as it tra¬ 
verses a particle detector at the LHC. In general gauge- 
mediated supersymmetry breaking (GGM) scenarios [4], 
the next-to-lightest supersymmetric particle (NLSP) de¬ 
cays into an SM particle and the LSP, which is a very 
light gravitino. The NLSP width is suppressed by the 
large supersymmetry-breaking scale, and may be such 
that its decay leads to the formation of a DV. Within 
split supersymmetry [5, 6], gluino ( g ) decay is sup¬ 
pressed by the high mass of the squarks. Long-lived 
gluinos then hadronize into heavy “A-hadrons” that may 
decay at a detectable distance from their production 
point. Additional scenarios with LLPs include hidden- 
valley [7], dark-sector gauge bosons [8], and stealth su¬ 
persymmetry [9]. Some of the models are disfavored [10] 
by the recent observation of a Higgs boson at mu ss 
125 GeV [11, 12]. 

This article presents the results of a search for DVs 
that arise from decays of long-lived, heavy particles, at 
radial distances of millimeters to tens of centimeters from 
the proton-proton IP in the ATLAS detector at the LHC. 
Two types of signatures are considered. In the dilepton 
signature, the vertex is formed from at least two lepton 
candidates (with “lepton” referring to an electron or a 


muon), with opposite electric charges. In the multitrack 
signature, the DV must contain at least five charged- 
particle tracks. This signature is divided into four dif¬ 
ferent final states, in which the DV must be accom¬ 
panied by a high-transverse-momentum (high-px) muon 
or electron candidate that originates from the DV, jets, 
or missing transverse momentum (A™ lss ). These sig¬ 
natures are labeled DV+muon, DV+electron, DV+jets, 
and DV+A™ ISS , respectively. In all signatures, at least 
one DV is required per event. In all cases, the expected 
background is much less than one event. 

The multitrack results improve on the previous ATLAS 
searches for this signature [13, 14] in several ways. The 
LHC center-of-mass energy is increased to 8 TeV, and 
the integrated luminosity is more than 4 times larger. 
While the previous search required only a high-px muon 
trigger, the current search also uses high-px electron, jets, 
or A™ lss triggers. Furthermore, the detector volume used 
for the search has been enlarged by more than a factor 
of 3. 

This is the first search for high-mass, displaced lep¬ 
ton pairs at ATLAS. A previous ATLAS search [15] con¬ 
sidered pairs of muons that were highly collimated due 
to the low mass of the decaying particle. ATLAS has 
also searched for long-lived particles that decay inside 
the hadronic calorimeter [16, 17], the inner detector or 
the muon spectrometer [18], or that traverse the entire 
detector [19]. 

Related searches have been performed at other exper¬ 
iments. The CMS Collaboration has searched for de¬ 
cays of a long-lived particle into a final state contain¬ 
ing two electrons, two muons [20, 21], an electron and 
a muon [22], or a quark-antiquark pair [23]. The LHCb 
Collaboration has searched for long-lived particles that 
decay into jet pairs [24]. The Belle Collaboration has 
searched for long-lived heavy neutrinos [25], and the 
BABAR Collaboration has searched for displaced ver¬ 
tices formed of two charged particles [26] . The DO Collab¬ 
oration has searched for displaced lepton pairs [27] and 
bb pairs [28], and the CDF Collaboration has searched 
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for long-lived particles decaying to Z bosons [29]. LLPs 
have also been searched for by the ALEPH Collaboration 
at LEP [30], 

This article is organized as follows. First the ATLAS 
detector and event samples used are described in Secs. II 
and III, respectively. The event reconstruction and ver¬ 
tex selection criteria are given in Sec. IV, while the signal 
efficiency is detailed in Sec. V. The background estima¬ 
tion is given in Sec. VI, with the systematic uncertainties 
on background and signal in Sec. VII. Finally, the search 
results are given in Sec. VIII, along with their interpre¬ 
tations in various supersymmetric scenarios. 


II. THE ATLAS DETECTOR 

The ATLAS experiment [31] [63] is a multipurpose de¬ 
tector at the LHC. The detector consists of several layers 
of subdetectors. From the IP outwards there is an in¬ 
ner tracking detector (ID), electromagnetic and hadronic 
calorimeters, and a muon spectrometer (MS). 

The ID is immersed in a 2 T axial magnetic field and 
extends from a radius of about 45 mm to 1100 mm and 
to |z| of about 3100 mm. It provides tracking and vertex 
information for charged particles within the pseudorapid¬ 
ity region \p\ < 2.5. At small radii, silicon pixel layers 
and stereo pairs of silicon microstrip detectors provide 
high-resolution position measurements. The pixel system 
consists of three barrel layers, and three forward disks on 
either side of the IP. The barrel pixel layers, which are 
positioned at radii of 50.5 mm, 88.5 mm, and 122.5 mm 
are of particular relevance to this work. The silicon mi¬ 
crostrip tracker (SCT) comprises four double layers in 
the barrel and nine forward disks on either side. The 
radial position of the innermost (outermost) SCT barrel 
layer is 30.3 cm (52.0 cm). A further tracking system, a 
transition-radiation tracker (TRT), is positioned at larger 
radii, with coverage up to | p\ = 2.0. This device has two 
hit thresholds, the higher of which is used to assist in 
the identification of electrons through the production of 
transition radiation within the TRT. 

The calorimeter provides coverage over the range |? 7 | < 
4.9. It consists of a lead/liquid-argon electromagnetic 
calorimeter, a hadronic calorimeter comprising a steel 
and scintillator-tile system in the barrel region and a 
liquid-argon system with copper and tungsten absorbers 
in the end caps. 

The MS provides muon identification and contributes 
to the muon momentum measurement. This device has 
a coverage in pseudorapidity of \r]\ < 2.7 and is a three- 
layer system of gas-filled tracking chambers. The pseu¬ 
dorapidity region \r]\ < 2.4 is additionally covered by 
separate trigger chambers, used by the hardware trigger 
for the first level of triggering (level-1). The MS is im¬ 
mersed in a magnetic field that is produced by a set of 
toroid magnets, one for the barrel and one each for the 
two end caps. 

Online event selection is performed with a three-level 


trigger system. It is composed of a hardware-based 
level-1 trigger that uses information from the MS trigger 
chambers and the calorimeters, followed by two software- 
based trigger levels. 


III. DATA AND SIMULATED EVENTS 

The data used in this analysis were collected in 2012 
at a pp center-of-mass energy of y/s = 8 TeV. Af¬ 
ter the application of detector and data-quality require¬ 
ments, the integrated luminosity of the data sample is 
20.3 fb -1 . The uncertainty on the integrated luminosity 
is ±2.8%. It is derived following the same methodology 
as that detailed in Ref. [32]. With respect to the ori¬ 
gin of the ATLAS coordinate system at the center of the 
detector, the mean position of the pp collision, averaged 
throughout the collected data sample is ( x ) = —0.3 mm, 
( y ) = 0.7 mm, (z) = —7.7 mm. The RMS spread of the 
2 distribution of the collisions is a z = 47.7 mm, and the 
spreads in the x and y directions are less than 0.1mm. 

Samples of simulated Monte Carlo (MC) events are 
used to study the reconstruction and trigger efficiency 
for signal events within RPV, split supersymmetry, and 
GGM scenarios. In each simulated event, two gluinos or 
two squarks are created in the pp collision. Both of these 
primary particles undergo decay chains described by the 
same set of effective operators. In the simulated GGM 
and RPV scenarios, the LLP is the lightest neutralino 
Xi- In the split-supersymmetry scenario, the LLP is the 
gluino. Diagrams representing the simulated processes 
are shown in Fig. 1. 

All samples are generated with the AUET2B ATLAS 
underlying-event tune [33] and the CTEQ6L1 parton dis¬ 
tribution function (PDF) set [34], Events are generated 
consistently with the position of the pp luminous region 
and weighted so as to yield the correct 2 distribution of 
the collisions. Each generated event is processed with the 
GEANT4-based [35] ATLAS detector simulation [36] and 
treated in the same way as the collision data. The sam¬ 
ples include a realistic modeling of the effects of multiple 
pp collisions per bunch crossing observed in the data, ob¬ 
tained by overlaying additional simulated pp events gen¬ 
erated using PYTHIA 8 [37], on top of the hard scatter¬ 
ing events, and reweighting events such that the distri¬ 
bution of the number of interactions per bunch crossing 
matches that in the data. 

In what follows, the notation P —> A [L —> F] denotes 
an MC sample in which a primary particle P produced 
in the pp collision decays into a long-lived particle L and 
additional particles denoted A. The decay of the LLP 
into final state F is enclosed in square brackets. Samples 
where the primary particle is long-lived are denoted with 
[L —>• F], In both cases, masses may be indicated with 
parentheses, as in [L(100 GeV) — > F], The symbol q 
indicates a u- or d-quark unless otherwise specified, and 
£ indicates an electron or a muon. Charge conjugation of 
fermions is to be understood where appropriate. 
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RPV samples of the type g —> qq[xi —> M!v] are 
produced with HERWIG++ 2.6.3 [38]. Decays of the 
neutralino only into light leptons, which may be e + e _ , 
p + , or e ± /i =F , take place due to the nonzero values of 
the RPV couplings Am and Am [2]. 

RPV samples of q —► q[xi —» iqq/vqq] are generated 
with PYTHIA 6.426.2 [39]. The Xi decay into two light 
quarks and an electron, muon, or neutrino, is governed 
by the nonzero RPV coupling A' n . Samples contain¬ 
ing heavy-flavor quarks, q —>• q[x{ —t iqb] (produced 
with A' 13 0) and q —> g[x? -A £cb] (corresponding 

A' 23 7 ^ 0) are a l so generated, in order to study the 
impact of long-lived charm and bottom hadrons on the 
efficiency of DV reconstruction. A g — > qq[xi —> tqq] 
sample is used to quantify the effect of prompt NLSP de¬ 
cays on the reconstruction efficiency, by comparing with 
the corresponding model with squark production. 

PYTHIA 6.426.2 is used to produce GGM samples 
denoted g —► qq[x ? —> GZ ], in which the NLSP Xi is a 
Higgsino-like neutralino. Both the leptonic and hadronic 
decays of the Z boson are considered. 

Within a split-supersymmetry scenario, 
PYTHIA 6.427 is used to simulate production and 
hadronization of a primary, long-lived gluino. Geant4 
simulates the propagation of the A-hadron through 
the detector [40], and PYTHIA decays the A-hadron 
into a stable neutralino plus two quarks (u, d, s, c or 
b), a gluon, or two top quarks. The resulting samples 
are denoted [g -A qqx% [g -A gx% or [g -A ttx% 
respectively. 

Signal cross sections are calculated to next-to-leading 
order in the strong coupling constant, adding the re¬ 
summation of soft gluon emission at next-to-leading- 
logarithmic accuracy (NLO+NLL) [41-45]. The nomi¬ 
nal cross section and its uncertainty are taken from an 
envelope of cross section predictions using different PDF 
sets and factorization and renormalization scales, as de¬ 
scribed in Ref. [46]. 

In addition to these signal samples, MC samples of 
QCD dijet events, Drell-Yan events, and cosmic-ray 
muons are used for estimating some systematic uncer¬ 
tainties and some of the smaller background rates, as well 
as for validation of aspects of the background-estimation 
methods. 


IV. EVENT RECONSTRUCTION AND 
SELECTION 



FIG. 1: Diagrams representing some of the processes un¬ 
der study, corresponding to the simulated event samples. In 
RPV scenarios, the long-lived neutralino may decay via the 
(a) Xijk or (b) A [j k couplings, (c) Long-lived neutralino de¬ 
cay in a GGM scenario, (d) Long-lived A-hadron decay in a 
split-supersymmetry scenario. The quarks and leptons shown 
may have different flavors. Filled circles indicate effective in¬ 
teractions. 


candidates, jets, or A™ lss . This selection is described in 
Sec. IV A. 

Off-line selection criteria for leptons, jets, and A™ lss 
(see Sec. IV B) are used to further filter events for off¬ 
line processing, as described in Sec. IV C. 

Events satisfying the filter requirements undergo a 
CPU-intensive process termed “retracking”, aimed at ef¬ 
ficient reconstruction of tracks with large impact param¬ 
eter (do) with respect to the transverse position of any 
primary vertex (PV) of particles formed from the pp col¬ 
lision. Retracking is described in Sec. IV D. 

At the final event-selection stage, the presence of a pp 
collision is ensured by first requiring the event to have a 
PV formed from at least five tracks and situated in the 
longitudinal range \z\ < 200mm, consistent with the IP. 

The final selection is based on the reconstruction of a 
multitrack DV or dilepton DV, described in Secs. IV E 
and IV F, respectively. 


The event-reconstruction and selection procedures are 
designed, based on MC and experience from previous 
analyses [13, 14], to strongly suppress background and ac¬ 
commodate robust background-estimation methods (de¬ 
scribed in Sec. VI), while efficiently accepting signal 
events over a broad range of LLP masses, lifetimes, and 
velocities. 

The initial event selection is performed with a com¬ 
bination of triggers that require the presence of lepton 


A. Trigger requirements 

Events must satisfy trigger requirements based on 
muon, electron, jets, or E™' ss criteria. 

Where muon triggers are used, a muon candidate is 
required by the trigger algorithm to be identified in the 
MS with transverse momentum px > 50 GeV. Its pseu¬ 
dorapidity must be in the MS-barrel region r; < 1.07, 
to reduce the trigger rate from fake muons due to beam 
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background in the end cap region. 

Photon triggers are used for channels requiring elec¬ 
tron candidates, since the ID track of a high-do electron 
may not be reconstructed. These require only a high- 
energy deposit in the electromagnetic calorimeter, and 
have no veto or selection based on ID tracks. Photon trig¬ 
gers provide significantly less background rejection than 
muon triggers. Therefore, the trigger used for final states 
involving electrons requires either a single photon candi¬ 
date with pt > 120 GeV or two photon candidates with 
Pt > 40 GeV each. 

The trigger requirement for the DV+£’™ lss search is 
pjraiss > gQ Q e y The DV+jets trigger requires four jets 
with pt > 80 GeV, five jets with pt > 55 GeV, or six 
jets with pt > 45 GeV. 


B. Off-line object definition 

The reconstruction and selection criteria for E™ lss , 
muon, electron, and jet candidates are described in what 
follows. These object definitions are used by the off¬ 
line filter (Sec. IV C) and the final analysis (Secs. IV E 
and IV F). 


1. Muon selection 

Muon candidates are required to be reconstructed in 
both the MS and the ID. The ID track associated with 
the muon candidate is required to have at least four SCT 
hits, but the number of required hits is reduced if the 
track crosses nonoperational sensors. Furthermore, the 
track must satisfy an \r ]\-dependent requirement on the 
number of TRT hits. No pixel hit requirement is applied 
to the muon-candidate track, which is different from the 
standard ATLAS muon-reconstruction algorithm [47]. 


2. Photon and electron selection 

Photon and electron candidates are identified with cri¬ 
teria based on the fraction of the candidate energy de¬ 
posited in the hadronic calorimeter and on the shape of 
the electromagnetic shower. In addition, electron candi¬ 
dates must be in the pseudorapidity region \r/\ < 2.47, 
and must satisfy requirements on the number of TRT 
hits associated with the ID track, the fraction of high- 
threshold TRT hits, and the pseudorapidity difference 
between the electron-candidate track and the associated 
calorimeter cluster. In contrast to the standard ATLAS 
electron-selection requirements [48], no requirement on 
the number of silicon hits is applied. 


3. Jet and i?™ lss selection 

Jet candidates are reconstructed using the anti -kt jet 
clustering algorithm [49, 50] with a radius parameter 
R — 0.6. The inputs to this algorithm are the ener¬ 
gies of clusters [51, 52] of calorimeter cells seeded by 
those with energy significantly above the measured noise. 
Jet momenta are constructed by performing a four-vector 
sum over these cell clusters, treating each cell as a four- 
momentum with zero mass. Jets are initially calibrated 
to the electromagnetic energy scale, which correctly mea¬ 
sures the energy deposited in the calorimeter by elec¬ 
tromagnetic showers [51]. Further jet-energy scale cor¬ 
rections are derived from MC simulation and data, and 
used to calibrate the energies of jets to the scale of their 
constituent particles [51]. Jets are required to satisfy 
|? 7 | < 4.5 after all corrections are applied. 

A special category of jets is termed “trackless” jets. 
These are reconstructed as above, except that the anti-fct 
radius parameter is R = 0.4, the jet pseudorapidity is in 
the range \p\ < 2.5, and the scalar sum of the transverse 
momenta of the tracks in the jet is required to satisfy 
]T) tr pr < 5 GeV. Trackless jets may arise from decays 
of LLPs that take place far from the PV, where track- 
reconstruction efficiency is low. 

The measurement of the missing transverse momentum 
E™ lss is based on the calibrated transverse momenta of 
all jet and lepton candidates, as well as all calorimeter 
energy clusters not associated with such objects [53, 54], 


C. Off-line-filter requirements 

Events are selected for retracking and subsequent off¬ 
line analysis based on off-line filters that require one of 
the following: 

• A muon candidate with pt > 50 GeV, an elec¬ 
tron candidate with p T > 110 GeV, or a photon 
candidate with pt > 130 GeV. Electron candi¬ 
dates, and muon candidates that are associated 
with an ID track at this stage, are required to have 
do > 1.5 mm. The sample selected by this criterion 
contains 8.5 x 10 6 events. 

• A pair of candidate electrons, photons, or an 
electron-photon pair, with pt thresholds between 
38 and 48 GeV per object, and electron impact pa¬ 
rameter satisfying do > 2.0 mm or do > 2.5 mm 
depending on the channel. This criterion selects 
2.4 x 10 6 events. 

• Either two 50 GeV trackless jets and E™ 1SS > 
100 GeV (selecting 1.9 x 10 4 events), or one 45 GeV 
trackless jet and between four and six jets passing 
the same pt thresholds as those applied in the trig¬ 
ger, listed in Sec. IV A (selecting 4.6 x 10 5 events). 
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D. Retracking 

In standard ATLAS tracking [31], several algorithms 
are used to reconstruct charged-particle tracks. In the 
silicon-seeded approach, combinations of hits in the pixel 
and SCT detectors are used to form initial track can¬ 
didates (seeds) that are then extended into the TRT. 
Another algorithm starts from track segments formed of 
TRT hits, and extrapolates back into the SCT, adding 
any silicon hits that are compatible with the recon¬ 
structed trajectory. Both of these methods place con¬ 
straints on the transverse and longitudinal impact pa¬ 
rameters of track candidates that result in a low effi¬ 
ciency for tracks originating from a DV, many of which 
have large do- 

To recover some of these lost tracks, the silicon-seeded 
tracking algorithm is rerun off-line, using only hits that 
are not associated with existing tracks, for the events 
that satisfy the trigger and filter requirements (Secs. IV A 
and IV C). This retracking procedure is performed with 
the looser requirement do < 300 mm and |;?o| < 1500 mm. 
Furthermore, retracking requires a track to have at least 
five detector hits that are not shared with other tracks, 
while the corresponding requirement in standard silicon- 
seeded tracking is at least six hits. To reduce the rate 
of false seed tracks, it is required that these additional 
tracks have pt > 1 GeV, while the standard-tracking re¬ 
quirement is pt > 400 MeV. 

The remainder of the analysis proceeds with both the 
standard-tracking tracks and retracking tracks. To real¬ 
ize the benefits of retracking for lepton candidates, the 
lepton-identification algorithms are rerun with the re¬ 
tracking tracks. 

E. Multitrack vertex reconstruction and final 
selection 

1. Multitrack vertex reconstruction 

Tracks used for DV reconstruction are required to sat¬ 
isfy pt > 1 GeV and to have at least two SCT hits, to 
ensure high track quality. The requirement do > 2 mm is 
also applied, rejecting at least 97% of tracks originating 
from a PV. The tracks are rejected if they have no TRT 
hits and fewer than two pixel hits, in order to remove 
fake tracks. 

The selected tracks are used to construct a mul¬ 
titrack DV by means of an algorithm based on the 
incompatibility-graph approach [55]. 

The algorithm starts by finding two-track seed vertices 
from all pairs of tracks. Seed vertices that have a ver¬ 
tex fit y 2 of less than 5.0 (for 1 degree of freedom) are 
retained. If the seed vertex is inside the innermost pixel 
layer, both tracks must have a hit in this layer. If the 
vertex is between the first and second (second and third) 
pixel layers, both tracks must have a hit in either the 
second or third pixel layer (third pixel layer or the SCT). 


A seed vertex is rejected if any of its tracks have hits at 
radial positions smaller than that of the vertex. The in¬ 
teresting case of a charged LLP is not precluded by this 
selection, as the track formed by the LLP itself fails the 
do > 2 mm requirement, and is therefore not included in 
the seed vertex. To ensure consistency between the po¬ 
sition of the seed vertex and the direction p of the three- 
momentum vector of the seed-vertex tracks, the require¬ 
ment d ■ p > —20 mm is applied, where d = rbv — rpv is 
referred to as the “distance vector” between the position 
of the DV and that of the first PV. The first PV is de¬ 
fined as the PV with the largest Y^Pti where the sum is 
over tracks associated with the PV. 

Multitrack vertices are formed from combinations of 
seed vertices in an iterative process, as follows. If a track 
is assigned to several vertices, the vertex DVi with re¬ 
spect to which it has the largest y 2 is identified. If this 
X 2 is larger than 6, the track is removed from DVi. Oth¬ 
erwise, the algorithm finds the vertex DV 2 that has the 
smallest value of D/ctd , where D is the distance between 
DVi and DV 2 , and od is the estimated uncertainty on 
D. If D/ao < 3, a single vertex is formed from all the 
tracks of both vertices. If this is not the case, the track 
is removed from DVi. This process continues until no 
track is associated with more than one vertex. Finally, 
vertices are combined and refitted if they are separated 
by less than 1 mm. No requirement is made on the total 
charge of the tracks forming a vertex. 


2. Vertex selection 

The reduced y 2 of the DV fit is required to be smaller 
than 5.0. The DV position must be within the fiducial 
region rov < 300 mm, |zdv| < 300 mm, where r dv an d 
2 dv are the radial and longitudinal DV positions with 
respect to the origin. To minimize background due to 
tracks or iginating from the PVs, the tra nsverse distance 
A xy = sj (x D v - a:pv) 2 + (yov ~ 2/pv) 2 between the DV 
and any of the PVs is required to be at least 4 mm. Here 
x and y are the transverse coordinates of a given vertex, 
with the subscripts PV and DV denoting the type of 
vertex. 

DVs that are situated within regions of dense detec¬ 
tor material are vetoed using a three-dimensional map of 
the detector within the fiducial region. The map is con¬ 
structed in an iterative process, beginning with geomet¬ 
rically simple detector elements that are fully accounted 
for in the MC simulation. Subsequently, detailed struc¬ 
tures, as well as the positioning and thickness of the sim¬ 
ple elements, are obtained from the spatial distribution 
of vertices obtained from the data, taking advantage of 
the known <fi periodicity of the detector to reduce statisti¬ 
cal uncertainties. The vertices used to construct the map 
are required to be formed from fewer than five tracks, in 
order to avoid the signal region defined below. The in¬ 
variant mass of these vertices, assuming massless tracks, 
must be greater than 50 MeV, to suppress vertices from 
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photon conversions, which have low spatial resolution due 
to the small opening angle between the electrons, as well 
as electron scattering. Vertices arising from decays of 
K° s mesons are removed with an invariant-mass crite¬ 
rion. The transverse-plane projection of the positions of 
vertices that occur inside the material regions is shown 
in Fig. 2. 


10 3 


10 2 


10 


1 


FIG. 2: Transverse-plane density (in arbitrary units) of ver¬ 
tices with fewer than five tracks in material regions that are 
excluded by the material veto in the region \z\ < 300 mm. The 
innermost circle corresponds to the beam pipe. This is sur¬ 
rounded by the three pixel layers. The octagonal shape and 
outermost circles are due to support structures separating the 
pixel and SCT detectors. 

As the final step in multitrack DV selection, the num¬ 
ber of tracks forming the DV is required to satisfy _/V tr > 
5, and the invariant mass TOdv of all the tracks in the 
vertex to be greater than 10 GeV. In calculating todv, 
each track is taken to have the mass of the charged pion. 
Candidate vertices that pass (fail) the mov > 10 GeV 
requirement are hereafter referred to as being high-mpy 
(Iow-todv) vertices. 

The typical position resolution of the DV in the mul¬ 
titrack signal MC samples is tens of microns for 7 ~dv and 
about 200 /mi for zpv near the IP. For vertices beyond 
the outermost pixel layer, which is located at a radius 
of 122.5 mm, the typical resolution is several hundred 
microns for both coordinates. 


3. DV+lepton selection 

In the DV+muon search, the muon candidate is re¬ 
quired to have triggered the event and have transverse 


momentum px > 55 GeV, which is well into the re¬ 
gion where the trigger efficiency is approximately in¬ 
dependent of the muon momentum. The muon candi¬ 
date is further required to be in the range |p| < 1.07 
and have transverse impact parameter do > 1.5 mm. A 
cosmic-ray muon traversing the entire ATLAS detector 
is reconstructed as two back-to-back muon candidates. 
To reject cosmic-ray background, events are discarded 
if they contain two muon candidates with AI? cosm j c = 
1/(7 r — A cj)) 2 — (771 + 772) 2 < 0.04, where 771 and 772 are 
the pseudorapidities of the two reconstructed muon can¬ 
didates and A <f> is their angular separation in the az¬ 
imuthal plane. This has a negligible impact on the signal 
efficiency. 

In the DV+electron search, the electron candidate is 
required to have triggered the event and satisfy px > 
125 GeV and do > 1.5 mm. 

To ensure that the lepton candidate is associated with 
the reconstructed DV, the distance of closest approach 
of the selected muon or electron candidate to the DV is 
required to be less than 0.5 mm. This requirement en¬ 
sures that the reconstructed DV gave rise to the muon or 
electron candidate that triggered the event, and so the 
selection efficiency for each LLP decay is independent of 
the rest of the event. This facilitates a straightforward 
calculation of the event-selection efficiency for scenarios 
with different numbers of LLPs. The aforementioned se¬ 
lections are collectively referred to as the vertex-selection 
criteria. Events containing one or more vertices satisfy¬ 
ing these criteria are accepted. 


4- DV+jets and DV+E™ ISS selection 

The DV+jets selection requires one of the following: 
four jets with px > 90 GeV; five jets with px > 65 GeV; 
or six jets with px > 55 GeV. All jets considered in 
these selection criteria are required to have I 77 I < 2 . 8 . 
DV+jet candidate events are discarded if they contain 
any candidate jet failing to satisfy quality criteria de¬ 
signed to suppress detector noise and noncollision back¬ 
grounds [56, 57]. This has a negligible effect on the sig¬ 
nal efficiency. In the DV+F1™ 1SS search, the requirement 
p^miss > igQ Q e y j g applied. For these selection criteria, 
the trigger efficiency is approximately independent of the 
E“ lss and the jet transverse momenta. 


F. Dilepton selection 

In the dilepton search, muon candidates are required 
to satisfy px > 10 GeV, |? 7 | < 2.5, and do > 2mm. For 
electron candidates, the requirements are px > 10 GeV 
and do > 2.5 mm. A lepton candidate is discarded if 
its ID track is in the pseudorapidity region Ipl < 0.02, 
where the background-estimation procedure is observed 
to be unreliable (see Sec. VI). 
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To avoid double counting of vertices, lepton candidates 
used to form a dilepton DV must not have the same ID 
track as another lepton candidate. If two muon candi¬ 
dates or two electron candidates do share an ID track, 
the candidate that has the lower transverse momentum 
is discarded. If muon and electron candidates share an 
ID track, the electron candidate is discarded. Cosmic-ray 
muons, even those that interact while traversing the de¬ 
tector, are rejected by requiring that all lepton-candidate 
pairs satisfy Ai? cosm i c > 0.04. 

A dilepton DV is formed from at least two opposite- 
charge tracks identified as two electrons, two muons, or 
an electron and a muon. Any number of additional tracks 
may be included in the vertex. At this stage, it is verified 
that the dilepton selection criteria applied at the trigger 
and filter level (see Secs. IV A and IV C) are satisfied by 
the two lepton candidates forming the DV. Finally, the 
dilepton DV is required to satisfy the DV selection cri¬ 
teria specified in Sec. IV E 2, except for the requirement 
on the number of tracks, which is N tr > 2. As in the 
DV+lepton case, the dilepton-DV selection relies only 
on the leptons in the DV, and is independent of the rest 
of the event. 


V. SIGNAL EFFICIENCY 

In the dilepton and DV+lepton searches, where the 
selection criteria rely only on the particles produced in 
the DV, the vertex-level efficiency cdv is defined to be the 
product of acceptance and efficiency for reconstructing 
one signal DV, produced in the given search model, with 
all the trigger, filter, and final selection criteria. The 
event-level efficiency e e v defined as the probability for 
an event containing two DVs to be identified with at 
least one DV satisfying all the selection criteria, is then 
obtained from the relation 


The efficiency for reconstructing a multitrack DV is re¬ 
duced when the LLP decays to charm or bottom hadrons, 
resulting in two or more nearby DVs. Each of these DVs 
has a high probability of failing to meet the N tr and todv 
criteria, resulting in low efficiency if these DVs are not 
merged. This happens less at large values of rov; where 
DVs are more readily merged due to the worse position 
resolution. 
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e ev — 2 Z?£dv — ^ 2 £dv5 (1) 

where B is the LLP branching fraction into the specific 
search channel. In the DV+jets and DV+.E™ 1SS searches, 
only the event-level efficiency is defined, since the selec¬ 
tion criteria involve the entire event. 

The efficiency for reconstructing a multitrack or dilep¬ 
ton DV with the above selection criteria depends strongly 
on the efficiencies for track reconstruction and track se¬ 
lection, which are affected by several factors: (1) The 
number of tracks originating from the DV and their total 
invariant mass increase with the LLP mass. (2) More 
tracks fail the minimal-do requirement for small r^v, or 
when the LLP is highly boosted. (3) The efficiency for re¬ 
constructing tracks decreases with increasing values of do- 
(4) When an LLP decays at a radius somewhat smaller 
than that of a pixel layer, many tracks share hits on that 
pixel layer, failing to meet the track-selection criteria. 
The resulting impact on efficiency can be seen in Fig. 3 
at radii around 45 mm, 80 mm, and 115 mm. 


(b) 

FIG. 3: Comparisons of the vertex-level efficiency as a func¬ 
tion of the vertex radial position rnv for different RPV sam¬ 
ples. The vertical gray lines show the position of the first, 
second, and third pixel layers, (a) For the q —> q[xi /r qq] 
(A 211 ) samples with rriq = 700 GeV, comparing two cases of 
different LLP masses: m^o = 494 GeV, and rm ~0 = 108 GeV. 

(b) For the q ->■ q[x° -t /iqq] (A' 211 ) and q ->• q\x° ->• fiqb] 
(A213) samples, indicated by the relevant nonzero RPV cou¬ 
plings. 

The vertex-level efficiency does not depend appreciably 
on whether the primary particle is a squark or a gluino. 
However, the nature of the primary particle determines 
the number of jets, and hence impacts the event-level 
efficiency in the DV+jets and DV+E“ 1SS channels. 

Examples of the impact of LLP boost, mass, and 
heavy-flavor decays on the vertex-level efficiency are 
shown in Fig. 3 for q —» q[xi -+ fJ-qq] samples. As an 





example of the benefits of retracking, it is worthwhile to 
note that without retracking, the vertex-level efficiency 
for the m(x ?) = 494 GeV, \' 2 n 0 sample shown in this 

figure is about 1% at tdv = 80 mm, and is negligible for 
larger radii. 

Events in each MC sample are generated with a fixed 
value of the LLP lifetime tmc- To obtain the vertex-level 
efficiency for a different lifetime r, each LLP is given a 
weight 

Wbv(t, t) = — exp (— -, (2) 

T \tmc tJ 

where t is the true proper decay time of the gener¬ 
ated LLP. The vertex-level efficiency is then the sum 
of weights for LLPs that satisfy all the criteria in the 
sample. The same procedure is applied when calculating 
the event-level efficiency, except that the entire event is 
weighted by 

Wevt(*l,* 2 ,T) = Wbv(tl,T)Wbv(*2,T), (3) 

where t\ and £2 are the true proper decay times of the two 
LLPs in the event. Examples of the resulting dependence 
of £dv and e ev on the average proper decay distance ct 
are shown in Fig. 4. For most models considered in this 
analysis, the peak efficiency is typically greater than 5%, 
and it occurs in the range 10 < ct 100 mm. 

VI. BACKGROUND ESTIMATION 

The expected number of background vertices is esti¬ 
mated from the collision data for each channel. Since 
the number of events satisfying the final selection crite¬ 
ria is very small, the general approach is to first obtain 
a high-statistical-precision assessment of the probability 
for background-vertex formation using a large data con¬ 
trol sample. That probability is then scaled by the size of 
the signal-candidate sample relative to that of the control 
sample. 

A. Multitrack-vertex background estimation 

Background vertices that are due to accidental spa¬ 
tial crossing of tracks in a jet, particle interactions with 
material, or heavy-flavor decays have low values of todv 
and/or N tI and thus fail the selection requirements. Such 
vertices may contribute to high-?TiDv> high-A tr back¬ 
ground vertices via two mechanisms. 

• The dominant source of backgrounds are Iow-todv 
vertices that are accidentally crossed by an unre¬ 
lated, high-p T track at large angle [0( 1 radian)] to 
the other tracks in the vertex. This is referred to 
as the accidental-crossing background. 

• A much smaller background contribution is due to 
merged vertices. In this case, two Iow-todv vertices 




FIG. 4: (a) The event-level efficiency as a function of 

ct for split-supersymmetry [g — > g/qqXi( 100 GeV)] samples 
with various gluino masses, reconstructed in the DV+E“ ISS 
channel, (b) The vertex-level efficiency for the RPV g —► 
qq[Xi —» egv] samples with combinations of gluino and neu- 
tralino masses, reconstructed in the eg dilepton channel. The 
total uncertainties on the efficiencies are shown as bands (see 
Sec. VII). 

are less than 1 mm apart, and thus may be com¬ 
bined by the vertex-reconstruction algorithm into 
a single vertex that satisfies the iV tr and Todv cri¬ 
teria. 

The expected background levels from the two sources 
are estimated from the data. In order for the background 
estimate to have high statistical precision, it is performed 
with a large sample containing all events that have un¬ 
dergone retracking. This includes the events selected for 
this search, as described in Sec. IV C, as well as events 
used for other ATLAS analyses. The sample is divided 
into three subsamples, referred to as the muon stream, 
the electron stream, and the jets+E™ lss stream, with the 
name indicating the type of trigger used to select the 
events. The background level is estimated separately in 
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each of these streams with the methods described below, 
and the results are used for the DV+muon, DV+electron, 
and DV+jets and DV+P™ lss signal regions, respectively. 

To obtain the final background estimate in the signal 
region, the background estimate in each stream is multi¬ 
plied by a final-selection scale factor F ‘ stream ; which is the 
fraction of events in the given stream that satisfy the final 
event-selection criteria, other than the DV selection crite¬ 
ria. The values of these fractions are 0.08%, 5.0%, 1.45%, 
and 0.04%, for the DV+muon, DV+electron, DV+jet, 
and DV+i?™ lss searches, respectively. This use of p stream 
assumes that the average number of vertices per event, 
i\dj v , is independent of the selection criteria. Based on 
the change in iVwhen the final selection criteria are 
applied, an upward bias correction of 60% is applied to 
the estimated background level in the DV+electron chan¬ 
nel (the correction is included in the p stream = 5.0% 
value quoted above). The other channels have negligi¬ 
ble bias. A 10% systematic uncertainty is estimated for 
all channels from the statistical uncertainty of the bias 
estimate. 


1. Background from accidental vertex-track crossings 

The accidental-crossing background is estimated sepa¬ 
rately in six radial regions, ordered from the inside out. 
Region 1 is inside the beam pipe. Regions 2, 3, and 4 
correspond to the volumes just before each of the three 
pixel layers. Regions 5 and 6 are outside the pixel lay¬ 
ers. Region 5 extends outwards to rov = 180 mm, where 
there is essentially no detector material, while Region 6 
covers the volume from 180 < rov < 300 mm. In each 
region, a study of the todv distribution of V tr -track ver¬ 
tices, where N tI = 3 through 6, leads to identification of 
two types of background vertices, as follows. 

The first type, which dominates the low-mov range, 
is due to accidental track crossings in Region 1, and 
particle-material interactions in the other regions. This 
contribution to the todv spectrum is referred to as 
collimated-tracks background, reflecting the typically 
small angle between the tracks. The Todv distribution 
P c N f(m DV ) for this contribution is modeled from the N tr - 
track vertices for which the average three-dimensional an¬ 
gle between every pair of tracks is less than 0.5. In Fig. 5, 
P 3 co Vdv) is seen to fully account for 3-track vertices 
with todv less than about 3 GeV. However, it does not 
account for vertices with higher masses, particularly the 
signal region, todv > 10 GeV. 

The high-mDV part of the mov distribution is 
dominated by the second contribution, referred to as 
“DV+track”. In this case, a ( N tI — l)-track vertex is 
crossed by an unrelated track at a large angle with re¬ 
spect to the momentum vector of the vertex tracks. To 
construct a model of the DV+track m dv distribution of 
V tr -track vertices, every (7V tr — l)-track vertex, referred 
to as an “acceptor” vertex, is paired with a “donated” 
track that is taken from a “donor” vertex in another 


event. This is done for N tI — 1 in the range 2-5, where ac¬ 
ceptor vertices with five tracks are required to have mass 
below 10 GeV, to avoid the signal region. 

The pairing of a vertex and a track is performed with 
the following procedure. The donor vertex must satisfy 
all the DV selection criteria, except that the requirement 
on its mass is not applied, and it may have as few as 
two tracks. To ensure that the donated track is able to 
accurately model the effects of a large-angle crossing, it 
is required that the donor vertex be from the same radial 
region as the acceptor vertex, and that there is a large 
angle between the direction of the donated track and the 
distance vector of the donor vertex. 

In all regions apart from Region 1, the momentum vec¬ 
tor of the donated track is then rotated, so that its az¬ 
imuthal and polar angles (A^donor and Addonor) with re¬ 
spect to the distance vector of the acceptor vertex match 
those that it originally had with respect to the donor ver¬ 
tex. This ensures that the contribution of the donated 
track to the acceptor vertex mass correctly reflects the 
accidental-crossing probability as a function of Arj^ onol 
and A+ionor ■ 

Then, the four-momentum of the acceptor vertex and 
the rotated four-momentum of the track are added, ob¬ 
taining the todv value of the A tr -track vertex that would 
have been formed from an accidental crossing of the ac¬ 
ceptor vertex and the rotated donated track. The re¬ 
sulting todv distribution for the N pa ; rs DV+track pairs 
found in each region is denoted hN tT (mr>v), such that 

/•OO 

/ h NtT (mnv)dmr>v = V pairs . (4) 

Jo 

Tracks from donor vertices in Region 1 are treated differ¬ 
ently, since they tend to have high pseudorapidity, which 
impacts their DV-crossing probability more than their 
Ac^donor and A^donor values. Therefore, a Region-1 track 
is not rotated before its four-momentum is added to that 
of the acceptor vertex. 

The high-mov distribution for A tr -track vertices is 
then modeled by 

-P/v t r( m Dv) = /^v t r( m Dv)> (5) 

where 

^ylOGeV 

fw GeV ^3 ( m DV ) dm DV 

is the scale factor that normalizes the model to the 
data, and N^ 0GoV is the number of 3-track vertices with 
todv > 10 GeV. The model-predicted number of N tI - 
track background vertices with Todv > 10 GeV for a 
given stream and region is given by 

/>00 

^stream = / (m DV )dm D V (7) 

J 10 GeV 
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The model describes the high-mov background dis¬ 
tribution in data well, as seen in Fig. 5 for jets+-Kj? lss - 
stream 3-track and 4-track vertices in Region 6. Also 
shown is the collimated-track contribution, which ac¬ 
counts for the Iow-todv part of the distribution. Using 
4-track vertices to validate Eq. (7), the prediction for 
each of the three streams and six regions is compared 
with the observed number of vertices. The comparison, 
summarized in Fig. 6, shows good agreement within the 
statistical precision. 




FIG. 5: The mass distribution for (a) 3-track and (b) 

4-track vertices (data points) from the jets+.Ex 1SB stream 
in Region 6, overlaid with the model //i 3 (mDv) of Eq. (5) 
(yellow-shaded histogram) at high mass. The lower panel 
of each plot shows the ratio of the data to this model. The 
model for the collimated-track contribution P 3 o11 (?tidv) (blue- 
shaded histogram), which is correlated with the low-mass data 
but not used for estimating the signal-region background, is 
also shown. 


The final numbers of expected background vertices, af¬ 
ter multiplying 7V^ t r ® am by the scale factor p stream ; are 
shown in Table I. 



FIG. 6: Summary of the number of observed (black open 
circles) 4-track, high-mDv vertices in each of the radial regions 
and filter-selection streams and the predicted number (red 
triangles), from Eq. (7). In Region 1, the prediction includes 
the contribution from merging of two 2-track vertices (see 
Sec. VIA 2). The error bars on the prediction are too small 
to be visible, and in some bins no events are observed. 


TABLE I: Estimated numbers of background vertices sat¬ 
isfying all of the multitrack signal selection criteria, which 
arise from a low-mass DV accidentally crossed by an unre¬ 
lated track. In each entry, the first uncertainty is statistical, 
and the second is systematic (see Sec. VII). 


Channel 

No. of background vertices (xlO 3 ) 

DV+jet 

DV+Ex iss 

DV+muon 

DV+electron 

410 ± 7 ± 60 

10.9 ±0.2 ±1.5 

1.5 ±0.1 ±0.2 

207 ± 9 ± 29 


2. Background due to merged vertices 

In the last step of DV reconstruction (see Sec. IV E 1), 
vertices are combined if they are separated by less than 
1 mm. To estimate the background arising from this pro¬ 
cedure, the distribution of the distance g? 2 dv between 
two 2-track or 3-track vertices is studied. Each of the 
selected vertices is required to satisfy the DV selection 
criteria of Sec. IV E 2 except the ?bdv and N tr require¬ 
ments, and their combined mass is required to be greater 
than 10 GeV. To obtain a sufficient number of vertices 
for studying the c^dv distribution, the distribution is re¬ 
constructed from a much larger sample of vertex pairs, 
where each vertex in the pair is found in a different event. 
This is referred to as the “model” sample. 
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To validate the c? 2 DV distribution of the model sample, 
it is compared to that of vertices that occur in the same 
event, referred to as the “same-event” sample. It is found 
that the z positions of vertices in the same event are 
correlated, since more vertices are formed in high-track- 
multiplicity regions corresponding to jets. This effect is 
absent in the model sample. As a result, the distribu¬ 
tions of the longitudinal distance between the vertices in 
the model and the same-event samples differ by up to 
30% at low values of ci 2 DV- To correct for this differ¬ 
ence, each vertex pair in the model sample is weighted so 
as to match the 2 component distribution of the same- 
event sample. After weighting, the model distribution 
of the three-dimensional distance d 2 DV agrees well with 
that of the same-event sample in the entire study range of 
^ 2 DV < 120 mm. This is demonstrated in Fig. 7 for pairs 
of 2-track vertices and for the case of a 2-track vertex 
paired with a 3-track vertex. 

The background level for the analysis requirement of 
N tI > 5 tracks is estimated from vertex pairs where one 
vertex has two tracks and the other has three tracks. The 
area under the model distribution in the range d 2 DV < 
1 mm yields a background prediction of 0.02±0.02 events 
in each of the DV+lepton channels, and 0.03±0.03 events 
in the DV+jets and DV+.E™ 1SS channels. After multi¬ 
plication by F stleam , this background is negligible rela¬ 
tive to the accidental-crossing background, described in 
Sec. VIA 1. Background from the merging of two 3-track 
vertices or a 2-track and a 4-track vertex is deemed much 
smaller still. 


B. Dilepton-vertex background estimation 
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Background DVs in the dilepton search may arise from 
two sources: 

• The dominant background is due to accidental spa¬ 
tial crossings of unrelated lepton candidates that 
happen to come close enough to satisfy the vertex- 
reconstruction criteria. 

• Minor backgrounds, due to tracks originating from 
the PV wrongly associated with a DV, decays of 
SM long-lived particles, or cosmic-ray muons. The 
levels of background from these sources are deter¬ 
mined to be negligible relative to the accidental¬ 
crossing background. 

1. Background from accidental lepton crossing 

The level of the accidental-crossing background is esti¬ 
mated by determining the crossing probability, defined as 
the probability for two unrelated lepton-candidate tracks 
to be spatially nearby and reconstructed as a vertex. 
Pairs of opposite-charge lepton candidates are formed, 
where each lepton candidate in a pair is from a differ¬ 
ent event and satisfies the lepton-selection criteria. The 


FIG. 7: The distribution of the distance d 2 DV between (a) 
two 2-track vertices and (b) a 2-track vertex and a 3-track ver¬ 
tex with a combined mass above 10 GeV for the jets+i?“ lss 
stream data (data points) and in the model sample, in which 
the two vertices are in different events (histogram). A con¬ 
servative 100% uncertainty on the model is shown in the 
data/model ratio plot. The inset shows the o^dv distance up 
to values of 120 mm. The merged-vertex background estimate 
is determined from the area under the model distribution in 
the range d 2 DV < 1mm in (b). 


momentum vector of one of the two lepton candidates, 
selected at random, is rotated through all azimuthal an¬ 
gles by a step 5<j). At each rotation step, the two lepton 
candidates are subjected to a vertex fit and the DV se¬ 
lection criteria. If the pair satisfies the selection criteria, 
it is assigned a weighted probability 8(f)/2ir. Averaging 
the weighted probabilities over all pairs gives the prob¬ 
ability for a lepton-candidate pair to accidentally form 
a vertex. The probability is observed to be independent 
of Scj) for 5(f> < 0.03. To obtain the final background 
estimate, this probability is multiplied by the number 
Nu of data events containing two opposite-charge lep¬ 
ton candidates that satisfy the lepton-selection criteria. 
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TABLE II: Estimated numbers of background vertices satis¬ 
fying all of the dilepton signal selection criteria, arising from 
random combinations of lepton candidates. In each entry, the 
first uncertainty is statistical, and the second is systematic 
(see Sec. VII). 


Channel 

No. of background vertices (xlO ,! ) 

e + e“ 

p + p,~ 

1.0 ±0.2 tgie 

2.4 ±0.9 

2.0 ±0.5 


This procedure yields the background predictions shown 
in Table II. Compared with these predictions, the back¬ 
ground level for a 3-track vertex, where at least two of 
the tracks are lepton candidates, is negligible. 

The validity of this method for estimating the num¬ 
ber of dilepton DVs, including the assumption that the 
tracks are uncorrelated, is verified in several ways. Using 
Z —► p + p~ and tt MC samples, the procedure is applied 
to vertices formed from two lepton candidates, a lepton 
candidate and another track, or two tracks that are not 
required to be lepton candidates. It is observed that the 
method correctly predicts the accidental-crossing back¬ 
ground to within about 10%. The background-estimation 
method is tested also on pairs of tracks in the data, ex¬ 
cluding pairs of lepton candidates, with a variety of se¬ 
lection criteria. The predicted and observed numbers of 
background vertices are again found to agree to within 
10% for all selection criteria. The method also repro¬ 
duces well the distributions of roov, cdv, £dv, d ■ p, 
and the azimuthal angle between the two lepton candi¬ 
dates, in both MC simulation and data. As an example, 
Fig. 8 shows the todv and tdv distributions observed 
for data vertices composed of two nonlepton tracks and 
the distributions predicted by pairing two tracks in dif¬ 
ferent events. Some differences between the model and 
the data are seen at certain radii (e.g. rov < 50 mm and 
250 < tdv < 270 mm), but these do not substantially 
affect the total number of DVs and are covered by the 
assigned systematic uncertainty (see Sec. VIIA 2). The 
prediction is accurate down to DV masses of 6 GeV, well 
below the DV selection criterion of 10 GeV. At smaller 
masses, contributions from other background sources be¬ 
come significant. 

This background-estimation method ignores the possi¬ 
bility of angular correlations between the leptons forming 
a background vertex. The associated systematic uncer¬ 
tainty is described in Sec. VIIA 2. 


2. Minor backgrounds 

Backgrounds from the following sources are found to be 
negligible relative to the accidental-crossing background, 
and are therefore neglected. 

A potential source of background is prompt production 
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FIG. 8: Distributions of the (a) vertex mass and (b) vertex 
position radius for vertices composed of two nonlepton tracks 
in the data sample (data points), and the predicted model dis¬ 
tribution obtained from vertices formed by combining tracks 
from two different data events (shaded histograms). The ra¬ 
tio of the data to the model distributions is shown below each 
plot. The gray bands indicate the statistical uncertainties for 
the predicted distributions. The inset shows the mass distri¬ 
bution in the low-mass region, elsewhere mov > 10 GeV is 
required. In (a), the highest bin shows the histogram over¬ 
flow. 


of hard lepton pairs, notably from Z —► £ + £~ decays. Re¬ 
quiring A xy < 4 mm and removing the Todv > 10 GeV 
requirement yields no dilepton-vertex candidates, so the 
data show no evidence for prompt background. There¬ 
fore, MC simulation is used to estimate the probability 
for leptons originating from Z —► £ + £~ decays to satisfy 
the minimum-do requirements, the probability for such 
leptons to satisfy the vertex requirements, and the prob¬ 
ability for a Z —> £ + £~ event to pass the analysis kine¬ 
matic requirements. Multiplying the product of these 
probabilities by the number of Z —> £ + £~ events pro¬ 
duced at ATLAS yields an estimate of 10 -5 Z —» p + 
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events and 10 -4 Z —> e + e _ events in the A xy < 4 mm 
sideband. Thus, the background from this source is neg¬ 
ligible. 

Background from cosmic-ray muons is studied with the 
AR cosm i c distribution of the two highest-px muon candi¬ 
dates in each event, which satisfy the selection criteria 
except the Ai? cosm i c > 0.04 requirement. The distribu¬ 
tion drops rapidly as Ai? cosm ; c increases, with the highest 
pair having Ai? cosm i c = 0.014. The pairs that also satisfy 
the DV selection criteria constitute less than 7% of this 
sample and have a similar Ai? cosm i c distribution, termi¬ 
nating at A.Rcosmic ~ 0.0045. Therefore, it is concluded 
that the rate for cosmic-ray background muons satisfying 
the AR cosm i c >0.04 requirement is several orders of mag¬ 
nitude below the accidental-crossing background. In the 
case of a partially reconstructed cosmic-ray muon cross¬ 
ing a reconstructed lepton candidate from a pp collision, 
the two tracks are uncorrelated and any contribution to 
the background is already accounted for in the results 
shown in Table II. 

Background from decays of known long-lived hadrons 
is studied from vertices in which only one track is required 
to be a lepton candidate. It is found to be negligible, due 
to the small probability for a hadron to be misidentified 
as a lepton candidate and the mass resolution of the de¬ 
tector. 


VII. SYSTEMATIC UNCERTAINTIES AND 
CORRECTIONS 

The dominant systematic uncertainties are those as¬ 
sociated with the efficiency for reconstructing displaced 
electrons and with the jet and A™ lss selection criteria. 
Since the background level is low, uncertainties on the 
background estimation have a minor effect on the results 
of the analysis. The methods for evaluation of the sys¬ 
tematic uncertainties are described in detail below. 


A. Background-estimation uncertainties 

1. Multitrack DV background uncertainties 

The choice of the ttidv >10 GeV mass range for de¬ 
termining the scale factor / (see Sec. VIAl), as well as 
differences between the Todv distribution of the vertices 
and that of the model, are a source of systematic uncer¬ 
tainty on the background prediction. To estimate this 
uncertainty, / is obtained in the modified mass ranges 
todv > 5 GeV and todv > 15 GeV. The resulting 10% 
change in the background prediction for DVs passing the 
final selection requirements is used as a systematic un¬ 
certainty. An additional uncertainty of 10% is estimated 
from the variation of F stream as the selection criteria are 
varied (see Sec. VIA 1). Compared with these uncertain¬ 
ties, the uncertainty on the much smaller merged-vertex 
background level is negligible. 


2. Dilepton background uncertainties 

The background-estimation procedure for the dilepton 
search (see Sec. VIB 1) normalizes the background to the 
number Nu of events containing two lepton candidates 
that could give rise to a DV that satisfies the selection 
criteria. Contrary to the underlying assumption of the 
background estimation, the two lepton candidates may 
be correlated, impacting their probability for forming a 
high-?n-DV vertex. To study the impact of such correla¬ 
tion, Nu is recalculated twice, placing requirements on 
the azimuthal angle between the two lepton candidates, 
A (f>u, of 0.5 < A 4>u < 7T and 0 < A(f>u < 7r —0.5. The re¬ 
sulting variation yields the relative uncertainty estimates 
on N u of and for the p + , e ± pff 1 and 

e + e _ channels, respectively. 

An uncertainty of 15% on the background prediction 
is estimated from the validation studies performed using 
MC simulation and data, described in Sec. VIBl. The 
resulting systematic uncertainties are shown in Table II. 


B. Signal-efficiency uncertainties and corrections 

1. Trigger efficiency 

The muon trigger efficiency is studied with a “tag-and- 
probe” method, in which the invariant-mass distribution 
of pairs of tracks is fitted to the sum of a Z —i p + 
peak and a background contribution. To reduce the back¬ 
ground, one of the muon candidates (the “tag”) is re¬ 
quired to be identified as a muon. The muon-trigger ef¬ 
ficiency is determined from the fraction of Z —p + p~ 
decays in which the other muon candidate (the “probe”) 
satisfies the trigger criteria. Based on the results of 
this study in data and MC simulation, a correction of 
Ae = —2.5% is applied to the MC-predicted trigger effi¬ 
ciency. A total uncertainty of a e = 1.7% is estimated by 
comparing the trigger efficiency as a function of the muon 
candidate pt in data and MC simulation, and by com¬ 
paring the results of the tag-and-probe method, applied 
to Drell-Yan MC, with MC generator-level information. 
Similar studies of the trigger selections used for the elec¬ 
tron channels lead to Ae = —1.5% and a e = 0.8% for 
the pt > 120 GeV photon trigger, and Ae = —0.5%, 
a e = 2.1% for the two-photon p T > 40 GeV trigger. The 
jets and E™ lss triggers are fully efficient after the off-line 
cuts. 


2. Off-line track-reconstruction efficiency 

The uncertainty associated with the reconstruction ef¬ 
ficiency for tracks that originate far from the IP is es¬ 
timated by comparing the decay radius distributions for 
K°s mesons in data and MC simulation. The comparison 
is carried out with the ratio 
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Pi(K° s ) 


N? ata (K%) 

ww s r 


( 8 ) 


where i = labels four radial regions between 

5 mm and 40 mm, and Nf a ' ta ‘^ MC {Kg) is the number 
of K% mesons in radial region i in data/MC simula¬ 
tion, obtained by fitting the two-track mass distributions. 
In the calculation of Nf /lc (Kg), the I\g candidates in 
MC are weighted so that their r; and p t distributions 
match those seen in the data. The ratio pi(Kg) is con¬ 
structed separately for |p| < 1 and r; > 1. The dif¬ 
ference Api(Kg) = Pi(Kg) — pi(Kg) quantifies the ra¬ 
dial dependence of the data-MC discrepancy. The dis¬ 
crepancy is largest in the outermost radial region, with 
Ap 4 (I\g) = —0.03 for |?y| < 1 and Ap 4 (Kg ) = —0.2 for 
1 77 1 > 1. The statistical uncertainties on pi(Kg) are neg¬ 
ligible compared to these discrepancies. 

To propagate this maximal discrepancy into a conser¬ 
vative uncertainty on the signal efficiency, DV daughter 
tracks are randomly removed from signal-MC vertices be¬ 
fore performing the vertex fit. The single-track removal 
probability is taken to be Ap 4 (Kg )/2 in each of the two 
eta\ regions. The resulting change in the DV efficiency 
is taken as the tracking-efficiency systematic uncertainty. 
This uncertainty is evaluated separately for each value of 
cr, and is generally around 1 %. 


3. Off-line lepton-identification efficiency 

The lepton-identification efficiency uncertainty is de¬ 
termined in ATLAS using Z —> ft decays, and is typically 
less than 1%. For this analysis, an additional uncertainty 
associated with identification of high-do leptons is evalu¬ 
ated. 

For muons, this is done by comparing a cosmic-ray 
muon simulation to cosmic-ray muon candidates in data. 
The events are required to pass the muon trigger and to 
have two muon candidates that fail the muon veto (see 
Sec. IVBl). The MC muons are weighted so that their 
77 and (f> distributions are in agreement with those of the 
data. Comparing the ratio of the muon candidate do 
distributions in data and in MC simulation yields a do- 
dependent efficiency correction that is between 1 % and 
2.5%, with an average value of 1.5%. The uncertainty as¬ 
sociated with this procedure is taken from the statistical 
uncertainty, and is 2 % on average. 

Unlike in the case of cosmic-ray muons, there is no 
easily identifiable, high-rate source of large-do electrons. 
Therefore, the performance of the simulation is validated 
by comparing the electron-identification efficiency e e (zo) 
as a function of the longitudinal impact parameter Zq of 
the electron candidate in data and MC simulation, mea¬ 
sured with the tag-and-probe method using Z —i e + e~ 
events. It is observed that e e (^o) is consistent in data and 
in MC simulation to better than 1 % for |zo| < 250 mm, 
beyond which there are too few events for an accurate 


measurement. Furthermore, the value of e e obtained 
with the tag-and-probe method is consistent to within 
1% with that determined from MC generator-level infor¬ 
mation. This data-MC agreement in e e (zo) is taken as 
an indication that the do dependence of the efficiency, 
e e (do), is also well described by simulation. In signal MC 
samples, the function e e (do) varies by about 10 % due to 
kinematic correlations between do and factors that affect 
the efficiency, such as the value of r DV and the boost of 
the LLP. To account for the possibility of an additional 
do dependence that may not be well simulated, a system¬ 
atic uncertainty of 10 % on the electron-identification effi¬ 
ciency is assigned. This conservative uncertainty weakens 
the upper limits by less than 1 %. 


4 . Jets and E™ lss reconstruction 

The impact on the signal efficiency of uncertainties in 
the jet-energy scale calibration and jet-energy resolution 
is evaluated following the methods described in Refs. [57] 
and Ref. [58], respectively. An additional uncertainty on 
the jet pt is evaluated for jets that originate from the de¬ 
cay of an LLP, by linearly parametrizing the pt mismea- 
surement in MC simulation as a function of tdv and zdv • 
The only significant dependence observed is a relative 
Pt mismeasurement of (4 ± 1) x 10 _5 (rDv/ mm), which 
is propagated to the jet-selection efficiency as a system¬ 
atic uncertainty. To account for possible mismodeling of 
trackless jets, an uncertainty is obtained by varying the 
requirement on Xur Pt for these jets. Systematic uncer¬ 
tainties in the E™' ss measurement are evaluated with the 
methods described in Refs. [53, 54] and propagated to the 
efficiency uncertainty. The impact of uncertainties in the 
simulation of initial-state radiation is estimated by vary¬ 
ing the pt distribution of the primary particles according 
to the distribution observed in MADGRAPH5 [59] sam¬ 
ples. The resulting efficiency uncertainty varies between 
2 % and 10 %. 


5. Multiple pp interactions 

The dependence of the reconstruction efficiency on the 
number of pp interactions per bunch crossing is studied 
by varying the average number of interactions per LHC 
bunch crossing in the simulation by 4%. This value re¬ 
flects uncertainties in the detector acceptance, trigger ef¬ 
ficiency, and modeling of additional pp interactions. The 
resulting relative uncertainty is typically of order 1 % or 
less. 


VIII. RESULTS 

Figure 9 shows the distribution of ttidv versus the 
number of associated lepton candidates in the selected 
data sample before the final selection requirements on 



15 


TABLE III: Model-independent 95% confidence-level upper 
limits on the visible cross section for new physics in each of 
our searches. 


Channel 

Upper limit on visible cross section [fb] 

DV+jet 

0.14 

DV+U™ iss 

0.15 

DV+muon 

0.15 

DV+electron 

0.15 

e + e“ 

0.14 

1 

3. 

+ 

3. 

0.14 

e ± ^ 

0.15 


these variables are applied. The distributions of todv 
versus the number of tracks in the vertex obtained for 
the multitrack-DV search are shown in Figs. 10 and 11. 
No events are seen in the signal region for any of the seven 
channels. In addition, no same-charge dilepton vertices 
are seen with ttidv > 10 GeV. The distributions ex¬ 
pected for some of the signal samples are also shown for 
comparison. 

Given the lack of a signal observation, 95% confidence- 
level upper limits on the total visible cross section for new 
physics are shown in Table III. 

Furthermore, for each of the physics scenarios con¬ 
sidered, 95% confidence-level upper limits on the signal 
yields and production cross sections are calculated for dif¬ 
ferent values of the proper decay distance ct of the LLP, 
and presented in the figures in this section. The lim¬ 
its are calculated using the CLs prescription [60] with 
the profile likelihood used as the test statistic, using the 
HistFitter [61] framework. Uncertainties on the signal ef¬ 
ficiency and background expectation are included as nui¬ 
sance parameters, and the CLs values are calculated by 
generating ensembles of pseudoexperiments correspond¬ 
ing to the background-only and signal-plus-background 
hypotheses. 

Since less than one background event is expected in all 
cases and no events are observed, the observed limits are 
very close to the expected limits. 

In the case of the dilepton and DV+lepton searches, 
where the trigger and reconstruction depend almost ex¬ 
clusively on the signal DV, upper limits on the number of 
vertices produced in 20.3 fb -1 of data are presented for 
each channel, accounting for the vertex-level efficiency 
at each value of ct. Figure 12 shows these number lim¬ 
its for the DV+lepton search signatures. The limits are 
given separately for different masses of the long-lived neu- 
tralino and the primary squark or gluino, as well as for 
different A% fc couplings, which give rise to light- or heavy- 
flavor quarks in the final state (see Sec V for a discus¬ 
sion of the heavy-flavor and mass impact on efficiency). 
The number limits for the dilepton search signatures are 
shown in Fig. 13 for the final states ee, /i/i, e/x, as well as 
for the combination of Z —> ee, Z —» /x/x and Z —> rr. 

In addition, limits on the production cross sections for 


events are presented for the different simulated scenarios. 

Figures 14 and 15 show cross section upper limits 
obtained with the dilepton-DV search. Upper limits 
are shown for gluino-pair production in the RPV sce¬ 
nario, with neutralino decays determined by the choice 
of nonzero RPV coupling Am or A 122 , as well as within 
the GGM scenario with leptonic decays of the Z boson. 
For example, the RPV scenario is excluded for gluino 
mass uig = 600 GeV, neutralino mass m^o = 400 GeV, 
and neutralino proper decay distance in the range 0.7 < 
cr < 3 x 10 5 mm. 

cross section upper limits obtained with the 
multitrack-DV search are shown in the remaining 
figures. These limits are calculated up to proper decay 
distances of cr = 1 m, to avoid inaccuracies associated 
with reweighting events to very high lifetimes when the 
efficiency depends on both LLPs in the event. 

Figure 16 shows the upper limits on the production 
cross section of two squarks in the RPV scenario, with 
different squark and neutralino masses, as well as differ¬ 
ent A 7 parameters governing the neutralino decay. These 
limits are obtained with the DV+jets search, which re¬ 
sults in tighter limits than the DV+lepton searches for 
this scenario. These results exclude a rriq = 1 TeV squark 
for m,~ 0 = 108 GeV and 2.5 < ct < 200 mm with either 

^ A. 1 

light- or heavy-quark neutralino decays. Figure 17 shows 
the cross section upper limits for gluino-pair production 
within the GGM scenario, using hadronic Z decays. The 
scenario is excluded, for instance, for nig = 1.1 TeV 
and to^o = 400 GeV in the proper decay distance range 
3 < ct < 500 mm. 

Figure 18 shows the upper limits on gluino-pair pro¬ 
duction cross section in the split-supersymmetry model. 
These limits are obtained from the results of the 
DV+J?“ 1SS and DV+jets searches. The sensitivity is 
greater for the cases with m^o = 100 GeV than for those 
with m^o = nig — 480 GeV, and the DV+if™ lss search 
performs better than DV+jets in these scenarios, exclud¬ 
ing nig < 1400 GeV in the range of proper decay lengths 
15 mm < ct < 300 mm. 

In Figs. 19 and 20, the region of gluino mass versus, 
proper decay distance that is excluded by these limits 
is shown. The limit for each point in parameter space 
is taken from the channel that is expected to yield the 
most stringent limit, which is DV+i?™ ISS for most points. 
For the region of parameter space where the sensitivity 
is greatest, 20 mm < ct < 250 mm and m^o = 100 GeV, 
gluino masses of rrig < 1500 GeV are excluded. This 
range of masses is comparable to or slightly larger than 
those excluded by prompt searches [62] , searches for long- 
lived R- hadrons stopped in the ATLAS calorimeter [17], 
or searches for stable, massive, charged particles [19]. 
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FIG. 9: The distribution of dilepton-vertex candidates in terms of the vertex mass versus the number of lepton candidates in 
the vertex, in the (a) fj, + p~, (b) eand (c) e + e _ search channels. The data distributions are shown with red ovals, the 
area of the oval being proportional to the logarithm of the number of vertex candidates in that bin. The gray squares show the 
<jj(600 GeV) — > qq[xi{50 GeV) — > ppv / epv / eev\ signal MC sample. The shape of the background rriDv distribution arises partly 
from the lepton-candidate pt requirements. The signal region defined by the two-lepton and tudv > 10 GeV requirements is 
indicated. 
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(b) 


FIG. 10: The distribution of (a) DV+muon and (b) 

DV+electron candidates in terms of the vertex mass versus 
the number of tracks in the vertex. The data distribution 
is shown with red ovals, the area of each oval being pro¬ 
portional to the logarithm of the number of vertex candi¬ 
dates in that bin. The gray squares show the <j(700 GeV) —>• 
q[Xi(494 GeV) —¥ Eqq] RPV signal MC sample. The signal 
region N tl > 5, m dv > 10 GeV is indicated. 


FIG. 11: The distribution of (a) DV+jets and (b) D V+E^ iss 
candidates in terms of the vertex mass versus the number of 
tracks in the vertex. The data distribution is shown with 
red ovals, the area of each oval being proportional to the 
logarithm of the number of vertex candidates in that bin. 
The gray squares show the g( 1.1 TeV) —> qq{\ ?(400 GeV) —> 
GZ] GGM signal MC sample in (a) and the [<j(1.4 TeV) —> 
Xi(100 G eV)qq/g] split-supersymmetry sample in (b). The 
signal region N tl > 5, m dv > 10 GeV is indicated. 



































18 



(a) 


(b) 


C/5 

> 

03 

O 

CD 

■o 

O T- 

O 

CD 

_Q 


10 7 


,6 


10' 


10 D 


10 


= 10 d 


CD 

Cl 

Cl 

3 


10* 



(c) 


(d) 


FIG. 12: RPV-scenario upper limits at 95% confidence level on the number of neutralinos in 20.3 fb -1 that decay into (a) nqq 
(with q indicating a u- or d-quark), (b) fiqb and /icb (indicated by the nonzero RPV couplings A 213 and A 223 ? respectively), (c) 
eqq , and (d) eqb and ecb (Ai 13 and Ai 2 3 , respectively). The upper limits account for the vertex-level efficiency for each value 
of the neutralino proper decay distance cr. The different curves show the results for different masses of the primary gluino or 
squark and of the long-lived neutralino, while the shaded bands indicate Tier variations in the expected limit. 
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(a) (b) 




(c) 


(d) 


FIG. 13: Upper limits at 95% confidence level on the number of neutralinos in 20.3 fb 1 that decay into (a) eev in the RPV 
model, (b) e^iv in the RPV model, (c) /.i/xi' in the RPV model, and (d) ZG in the GGM model. The upper limits account for 
the vertex-level efficiency for each value of the neutralino proper decay distance cr. The different curves show the results for 
different masses of the primary gluino and of the long-lived neutralino, while the shaded bands indicate ±1(7 variations in the 
expected limit. In some cases limits are terminated for cr < 1 mm due to limited statistical precision. 








20 





cx [mm] 

FIG. 15: The 95% confidence-level upper limits, obtained 
from the dilepton search, on the production cross section for a 
pair of gluinos of mass 1.1 TeV that decay into two quarks and 
a long-lived neutralino in the GGM scenario for two values of 
the neutralino mass.. For further details see Fig. 14. 


(b) 

FIG. 14: The 95% confidence-level upper limits, obtained 
from the dilepton search, on the production cross section for 
a pair of gluinos of different masses that decay into two quarks 
and a long-lived neutralino in different models: (a) the RPV 
scenario with a pure Am coupling, (b) the RPV scenario 
with a pure A 122 coupling. All relevant final-state lepton- 
flavor combinations are used. The shaded bands around the 
observed limits indicate ±1 a variations in the expected limit, 
while the horizontal bands show the theoretical cross sections 
and their uncertainties. In some cases limits are terminated 
for cr < 1 mm due to limited statistical precision. 
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(a) 


(b) 


FIG. 16: The 95% confidence-level upper limits, obtained from the DV+jets search, on the production cross section for a 
pair of squarks in the RPV scenario, with the neutralino decaying to a lepton and two quarks, according to the nonzero X' 
couplings indicated in each case. The squark and neutralino masses are (a) 700 GeV and 494 GeV or (b) 1 TeV and 108 GeV, 
respectively, representing the results for different neutralino masses and boosts. For further details see Fig. 14. 
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FIG. 17: The 95% confidence-level upper limits, obtained from the (a) DV+i?™ lss and (b) DV+jets searches, on the production 
cross section for a pair of gluinos of mass 1.1 TeV that decay into two quarks and a long-lived neutralino in the GGM scenario. 
For further details see Fig. 14. 
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(a) (b) 




(e) (f) 

FIG. 18: The 95% confidence-level upper limits, obtained from the (a, c, e) DV+-E™ 1SS and (b, d, f) DV+jets searches, on the 
cross section for gluino pair production in the split-supersymmetry model, with the gluino decaying to a neutralino plus either 
(a, b) a gluon or a light-quark pair or (c, d, e, f) a pair of top quarks. The mass of the neutralino is 100 GeV in (a, b, c, d) 
and is 480 GeV smaller than the gluino mass in (e, f). For further details see Fig. 14. 
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FIG. 19: The 95% confidence-level excluded regions fie below 
the curves shown in the mass-versus-cr plane for the split- 
supersymmetry samples, with the gluino decaying into a gluon 
or light quarks, plus a 100 GeV neutralino. The shaded bands 
indicate ±1 <t variations in the expected limit, while the dotted 
lines indicate the effect of varying the production cross section 
by l standard deviation. The expected and observed limits 
are identical. 




FIG. 20: The 95% confidence-level excluded regions lie below 
the curves shown in the mass-versus-cr plane for the split- 
supersymmetry samples, with the gluino decaying into (a) a 
top quark pair and a 100 GeV neutralino, or (b) a top quark 
pair and a neutralino with a mass that is 480 GeV smaller 
than the gluino mass. The shaded bands indicate ±lcr varia¬ 
tions in the expected limit, while the dotted lines indicate the 
effect of varying the production cross section by 1 standard 
deviation. The expected and observed limits are identical. 
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IX. SUMMARY AND CONCLUSIONS 

This article reports on a search for long-lived parti¬ 
cles decaying into two leptons or five or more charged 
particles. In the latter case, events are selected using as¬ 
sociated lepton candidates, jets or missing transverse mo¬ 
mentum. The main signature of the search is a displaced 
vertex with an invariant mass greater than 10 GeV. This 
signature corresponds to a wide variety of new-physics 
models, many of which have not been searched for pre¬ 
viously. The search uses the full data sample of pp colli¬ 
sions collected by the ATLAS detector at the LHC with 
a center-of-mass energy of y/s = 8 TeV and an integrated 
luminosity of 20.3 fb -1 . Less than one background event 
is expected in each of the channels, and no events are 
observed. Upper limits are provided on the number of 
long-lived particle decays in the data sample and on the 
cross section for production of particles that give rise 
to the search signatures in a variety of supersymmetric 
models. These upper limits exclude significant regions of 
the parameter space of new-physics models with particle 
masses within reach of the LHC. 
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